图像直方图概述、计算与绘制

文章目录

  • 1.图像直方图概述
  • 2.直方图的计算与绘制
  • 2.1 计算直方图:calcHist() 函数
  • 2.2 寻找最值:minMaxLoc() 函数
  • 2.3 示例程序:绘制 H-S 直方图
  • 2.4 示例程序:计算并绘制图像一维直方图
  • 2.5 示例程序:绘制 RGB 三色直方图

1.图像直方图概述


  直方图广泛运用于很多计算机视觉运用当中,通过标记帧与帧之间显著的边缘和颜色的统计变化,来检测视频中场景的变化。在每个兴趣点设置一个有相近特征的直方图所构成 “标签”,用以确定图像中的兴趣点。边缘、色彩、角度等直方图构成了可以被传递给目标识别分类器的一个通用特征类型。色彩和边缘的直方图序列还可以用来识别网络视频是否被复制。
  其实,简单来说,直方图就是对数据进行统计的一种方法,并且将统计值组织到一系列实现定义好的 bin 当中。其中, bin 为直方图中经常用到的一个概念,可以译为 “直条” 或 “组距”,其数值是从数据中计算出的特征统计量,这些数据可以是诸如梯度、方向、色彩或任何其他特征。且无论如何,直方图获得的是数据分布的统计图。通常直方图的维数要低于原始数据。
  图像直方图(Image Histogram)是用以表示数字图像中亮度分布的直方图,标绘了图像中每个亮度值的像素数。这种直方图中,横坐标的左侧为纯黑、较暗的区域,而右侧为较亮、纯白的区域。因此一张较暗图片的直方图中的数据多集中于左侧和中间部分,而整体明亮、只有少量阴影的图像则相反。CV 领域常借助图像直方图来实现图像的二值化。
  直方图的意义如下:
  ● 直方图是图像中像素强度分布的图形表达方式。
  ● 它统计了每一个强度值所具有的像素个数。

  直方图是对数据的统计集合,并将统计结果分布于一系列预定义的 bins 中。这里的数据不仅仅指的是灰度值,且统计数据可能是任何有效描述图像的特征。
  假设有一个矩阵包含一张图像的信息(灰度值 0 - 255),既然已知数字的范围包含 256 个值,于是可以按一定规律将这个范围分割成子区域(也就是 bins)。如:

  然后再统计每一个 bin(i) 的像素数目。采用这一方法来统计上面的数字矩阵,可以得到下图(其中 x 轴表示 bin,y 轴表示各个 bin 中的像素个数):

直方图的一些术语和细节:
  ● dims:需要统计的特征数目。在上例中,dims = 1 ,因为仅仅统计了灰度值(灰度图像)。
  ● bins:每个特征空间子区段的数目,可译为 “直条” 或 “组距”,在上例中, bins = 16。
  ● range:每个特征空间的取值范围。在上例中,range = [0, 255]。

2.直方图的计算与绘制
  直方图的计算在 OpenCV 中可使用 calcHist() 函数,而计算完成之后,可以采用 OpenCV 中的绘图函数,如绘制矩形的 rectangle() 函数,绘制线段的 line() 来完成。2.1 计算直方图:calcHist() 函数
  calcHist() 函数用于计算一个或多个阵列的直方图。
 

void calcHist(
	const Mat* images,	//输入的数组或数据集
	int nimages,		//输入数组的个数
	const int* channels,	//需要统计的通道(dim)索引	
	InputArray mask,	//可选的操作掩码,用于标记出统计直方图的数组元素数据
	OutputArray hist,	//输出的目标直方图	
	int dims,			//需要计算的直方图的维数
	const int* histSize,	//存放每个直方图尺寸的数组
	const float** ranges,	//每一维数值的取值范围
	bool uniform = true,	//指示直方图是否均匀的标识符
	bool accumulate = false		//累计标识符,主要是允许多从多个阵列中计算单个直方图,
								//或者用于在特定的时间更新直方图。
)

  ● 第一个参数:const Mat* 类型的 iamges,输入的数组(或数据集),它们需为相同的深度(CV_8U 或 CV_32F)和相同的尺寸。
  ● 第二个参数:int 类型的 nimages,输入数组的个数,也就是第一个参数中存放了多少张 “图像”,有几个原数组。
  ● 第三个参数:const int* 类型的 channels,需要统计的通道(dim)索引。第一个数组通道从 0 到 images[0].channels() - 1,而第二个数组通道从 images[0].channels() 计算到 images[0].channels() + images[1].channels() - 1。
  ● 第四个参数:InputArray 类型的 mask,可选的操作掩码。如果此掩码不为空,那么它必须为 8 位,并且与 images[i] 有同样的大小和尺寸。这里的非零掩码元素用于标记出统计直方图的数组元素数据。
  ● 第五个参数:OutputArray 类型的 hist,输出的目标直方图,一个二维数组。
  ● 第六个参数:int 类型的 dims,需要计算的直方图的维数,必须是正数,且不大于 CV_MAX_DIMS(在 OpenCV3 中等于 32)。
  ● 第七个参数:const int* 类型的 histSize,存放每个维度的直方图尺寸的数组。
  ● 第八个参数:const float** 类型的 ranges,表示每一个维度数组(第六个参数 dims)的每一维的边界阵列,可以理解为每一维数值的取值范围。
  ● 第九个参数:bool 类型的 uniform,指示直方图是否均匀的标识符,有默认值 true。
  ● 第十个参数:bool 类型的 accumulate,累计标识符,有默认值 false。若其为 true,直方图在配置阶段不会被清零。此功能主要是允许多从多个阵列中计算单个直方图,或者用于在特定的时间更新直方图。

原文查看转载文章

https://blog.csdn.net/qq_38701868/article/details/89215881#1_2

第一课:绘制直方图 import cv2 as cv import numpy as np from matplotlib import pyplot as plt def plot_demo(image): plt.hist(image.ravel(),256,[0,256]) plt.show() def image_hist(image): color = ('blue','green','red') for i,color in enumerate(color): #计算直方图 hist = cv.calcHist(image,[i],None,[256],[0,256]) #画出直方图 plt.plot(hist,color=color) plt.xlim(0,256) plt.show() src = cv.imread("E:/opencv/picture/Greatwall.jpg") cv.imshow("initial_window",src) plot_demo(src) image_hist(src) cv.waitKey(0) cv.destroyAllWindows() 效果图: 7-1源图像 7-2 plot_demo 7-3 image_hist 分析: 1. 什么叫直方图? 答:简单来说就是图像每个像素值的个数统计。比如说一副灰度图像素值为0的有多少个,像素值为1的有多少个….直方图是一种分析图片的手段。 在计算直方图之前,有几个术语需要先了解一下: dims:要计算的通道数,对于灰度图dims = 1,对于rgb图像dims= 3 range:要计算的像素值范围,一般为【0,256】(不包括256) bins:子区间数目,如果我们统计0~255每个像素值,bins =256;如果划分区间,比如【0,15】,【16,31】,….【240,255】这样的16个区间,bins = 16 2. 需要安装matplotlib pip install matplotlib 3. plot_demo(matplotlib自带的计算绘制直方图) def plot_demo(image): plt.hist(image.ravel(),256,[0,256]) plt.show() 1) image.ravel()函数的功能是将多维数组降为一维数组 2)matplotlib.pyplot.hist函数主要是计算直方图 hist函数原型:hist(x, bins=None, range=None,….) x参数表示是一个数组或一个序列,是指定每个bin(箱子)分布的数据 bins参数表示指定bin(箱子)的个数,也就是总共有几条条状图 range参数表示箱子的下限和上限。即横坐标显示的范围,范围之外的将被舍弃。 range这个参数一般都是[0,256] 3)plt.show()显示直方图 4. image_hist def image_hist(image): #注意这个是tuple()哦,且里面的是字符‘’而不是字符串 color = ('blue','green','red') for i,color in enumerate(color): hist =cv.calcHist(image,[i],None,[256],[0,256]) plt.plot(hist,color=color) plt.xlim(0,256) plt.show() 1) enumerate() 函数用于将一个可遍历的数据对象(如列表、元组或字符串)组合为一个索引序列,同时列出数据下标和数据,一般用在 for 循环当。 语法: enumerate(sequence, [start=0]) 参数: • sequence -- 一个序列、迭代器或其他支持迭代对象。 • start -- 下标起始位置。 返回值 返回 enumerate(枚举) 对象。 举个例子: 对于普通的for循环: >>>i = 0 >>> seq = ['one', 'two', 'three'] >>> for element in seq: print i, seq[i] i +=1 0 one 1 two 2 three 对于for循环使用enumerate: >>>seq = ['one', 'two', 'three'] >>> for i, element in enumerate(seq): print i, element 0 one 1 two 2 three 2)cv2.calcHist的原型为:calcHist(images, channels, mask, histSize, ranges[, hist[, accumulate]]) -> hist images参数表示输入图像,传入时应该用括号[ ]括起来 channels参数表示传入图像的通道,如果是灰度图像,那就不用说了,只有一个通道,值为0,如果是彩色图像(有3个通道),那么值为0,1,2,选择一个,对应着BGR各个通道。这个值也得用[ ]传入。 mask参数表示掩膜图像。如果统计整幅图,那么为None。主要是如果要统计部分图的直方图,就得构造相应的掩膜来计算。 histSize参数表示灰度级的个数,需要括号,比如[256] ranges参数表示像素值的范围,通常[0,256]。此外,假如channels为[0,1],ranges为[0,256,0,180],则代表0通道范围是0-256,1通道范围0-180。 3)plt.plot与plt.xlim函数 plt.plot(hist,color=color) plt.xlim(0,256) plt.plot函数:绘制函数图像 参数一:列表或者数组 参数二:控制颜色参数 plt.xlim函数:设置坐标轴刻度的取值范围 第二课 直方图均衡化与比较 均衡化(全局与部分) opencv的直方图均衡化都是基于灰度图像的 作用:直方图均衡化可以调整图像的对比度,使图像更加清晰,对比度增强。是图像增强的一个手段。 直方图均衡化:如果一副图像的像素占有很多的灰度级而且分布均匀,那么这样的图像往往有高对比度和多变的灰度色调。直方图均衡化就是一种能仅靠输入图像直方图信息自动达到这种效果的变换函数。它的基本思想是对图像像素个数多的灰度级进行展宽,而对图像像素个数少的灰度进行压缩,从而扩展像元取值的动态范围,提高了对比度和灰度色调的变化,使图像更加清晰。 源码: import cv2 as cv import numpy as np from matplotlib import pyplot as plt def hist_demo(image): hist = plt.hist(image.ravel(),256,[0,256]) plt.show() def Image_hist(image): color = ("b","g","r") for i,color in enumerate(color): #计算直方图 hist = cv.calcHist(image,[i],None,[256],[0,256]) plt.xlim(0,256) plt.plot(hist,color) plt.show() #opencv的直方图均衡化都是基于灰度图像的 #直方图均衡化可以自动的调整图像的对比度,使图像更加清晰化,对比度增强,是图像增强的一个手段 def equalHist_demo(image): gray = cv.cvtColor(image,cv.COLOR_BGR2GRAY) cv.imshow("gray",gray) dst = cv.equalizeHist(gray) cv.imshow("dst_win1",dst) #可以明显的看出dst_win1的对比度要比gray的对比度要高好多。 def clahe_demo(image): gray = cv.cvtColor(image, cv.COLOR_BGR2GRAY) cv.imshow("gray", gray) clahe = cv.createCLAHE(clipLimit=5.0,tileGridSize=(8,8)) dst = clahe.apply(gray) cv.imshow("dst_win2", dst) #之前学习的调节对比度以及锐化图像的方法 def contrast_demo(image,c,gamm): #blank = np.zeros(image.shape,np.uint8) #dst = cv.addWeighted(image,c,blank,0,gamm) kernel = np.array([[0,-1,0],[-1,5,-1],[0,-1,0]]) dst = cv.filter2D(image,-1,kernel) cv.imshow("dst_win",dst) src = cv.imread("E:/opencv/picture/lena.jpg") cv.imshow("inital_win",src) #hist_demo(src) #Image_hist(src) #contrast_demo(src,2,2) equalHist_demo(src) cv.waitKey(0) cv.destroyAllWindows() 得出结果如下: 原图 灰度图 图像全均衡化 自适应直方图均衡化(均衡参数可调推荐) 分析: 1)#opencv的直方图均衡化都是基于灰度图像的 #直方图均衡化可以自动的调整图像的对比度,使图像更加清晰化,对比度增强,是图像增强的一个手段 #全局直方图均衡化 def equalHist_demo(image): gray = cv.cvtColor(image,cv.COLOR_BGR2GRAY) cv.imshow("gray",gray) dst = cv.equalizeHist(gray) cv.imshow("dst_win1",dst) #可以明显的看出dst_win1的对比度要比gray的对比度要高好多。 api: cv.equalizeHist(src) 对输入的源灰度图像进行直方图均衡化处理,从而提高图像质量。如果源图像是彩色图像,得需先转化成灰度图像。 2) 自适应(局部)直方图均衡化(可调参数) .全局直方图均衡化可能得到是一种全局意义上的均衡化,但是有的时候这种操作并不是很好,会把某些不该调整的部分给调整了。Opencv还有一种直方图均衡化,它是一种局部直方图均衡化,也就是是说把整个图像分成许多小块(比如按10*10作为一个小块),那么对每个小块进行均衡化。 注:全图的直方图均衡化会导致对比度过度增强,所以在一些情况下应使用局部直方图均衡化; #局部直方图均衡化 def clahe_demo(image): gray = cv.cvtColor(image, cv.COLOR_BGR2GRAY) cv.imshow("gray", gray) clahe = cv
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值